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Abstract 

Certain invasive plants may rely on interference mechanisms (allelopathy, e.g.) 
to gain competitive superiority over native species. But expending resources on 
interference presumably exacts a cost in another life- history trait, so that the 
significance of interference competition for invasion ecology remains uncertain. 
We model ecological invasion when combined effects of preemptive and inter- 
ference competition govern interactions at the neighborhood scale. We consider 
three cases. Under "novel weapons," only the initially rare invader exercises 
interference. For "resistance zones" only the resident species interferes, and fi- 
nally we take both species as interference competitors. Interference increases 
the other species' mortality, opening space for colonization. However, a species 
exercising greater interference has reduced propagation, which can hinder its 
colonization of open sites. Interference never enhances a rare invader's growth 
in the homogeneously mixing approximation to our model. But interference 
can significantly increase an invader's competitiveness, and its growth when 
rare, if interactions are structured spatially. That is, interference can increase 
an invader's success when colonization of open sites depends on local, rather 
than global, species densities. In contrast, interference enhances the common, 
resident species' resistance to invasion independently of spatial structure, un- 
less the propagation-cost is too great. Increases in background mortality (i.e., 
mortality not due to interference) always reduce the effectiveness of interference 
competition. 
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1. Introduction 



Both lateral and vertical interactions can affect the likelihoo d that an inva- 

2003). But for 



sive species advances when rare (jLevine et al.l . I2004t iGoing et al 



many plants, competitive asymmetry between invader and re sident species gov- 

erns both the outcome and the timescale of eco l ogical i nvasion ( Lavergne and Molofskv 
2004 IVilaand Wchral, 12004 lO'Mallev etHI . l2006at iMacDougall et al.l l2009l )" 
Given invader-resident competition, eco l ogical superiority may depend on mor e 
than one mechanism (jCase and Gilpinl . Il974i iRidenour and Callawavl 1200 11 ) . 
Our analysis addresses the combined impact of preemptive and interference com- 
petition on invasion dynamics when bio tic interactions are structured spatially . 
We a s sume that int e rferen ce has a cost ( Adams et al. . Il979t iBazzaz and Graca 
1997t lAmarasekarel 120021 ): an increasing level of interspecific interference re- 
quires a reduction in propagation rate, diminishing that species' capacity to 
colonize unoccupied sites. 

In many plant communities the primary mode of competiti on is site preemp- 
tion (.Schoener . 1983; Bergelson. 1990; Crawlev et al. , 1999: Y urkonis and Meinerd . 
20041 ): i.e., species interact through colonization of empty sites. Superior pre- 
emptive competitors have higher propagat ion rates or lower mortality rates 
( Korniss and Caracol liool lO'Mallev et all . i2006b,: AUstadt et al.. . ,20071) . The 
former increases colonization of open sites, and the latter decreases a com- 
petitor's opportunities for colonization. Preemptive competitors have the same 
niche in a spatially homoge neous environment ( Amarasekare, 2003) . Ordinarily 



this precludes coexistence (jShurin et al.l . 12004 ; 



Allstadt et alH2009[). sin ce self 



Case and Gilph] ()l974h 



regulation does not exceed interspecific competition, 
suggest that this niche similarity might favor evolution of interference mecha- 
nisms. 

An interference competitor inhibits another species' access to a critical re- 
source, often by harming individuals of the other species. Examples i nclude in- 



terspecific territoriality in animal s and c hemical competition in plants (jCase et al 



1994 ICallawav and Ascheghouj. l2000t ). Exotic inva ders may suppress native 
species' densities through interference competition dP'Antoiiio and Vitousek , 



I992I : Callawav and Ridenour , 2004 ; Cappuccino and ArnasonL 2006[ ). The "novel 
weapons" hypothesis proposes tha t some in vasive plants release chemicals that 
inhibit growth of native species (jCallawav a nd Aschcghoug, 2000). Interest- 
ingly, allelopathic interference may act directly on individuals of the resident 
competitor, or may act indirectly through toxic effects on native species' mi- 
crobial mutuali s ts, pa rticularly mycorrhizal fungi ( Wolfe and Klironomosl 2005 : 



Callawav et al. . 20081 ). So, under the novel weapons hypothesis, invaders attain 



competitive superiority because their phytochemicals present novel challenges to 
native species. Reasonably, in other communities, exotic species likely en c ounte r 
interference competition from natives; see comments in I Von Holle et al. I (l2003h . 
Our models address the three general cases where interference can affect the 
outcome of resident-invader competition. We associate novel weapons with in- 
terference by the invader only. We refer to "resistance zones" when the resident 
species, but not the invader, exerts interference competition. And, of course. 
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invad ers and residents may each compete via interference (jCase and Gilpin . 
1974r ). If neither species exhibits interference, our model leaves preemption as 
the sole competitive mechanism. In this case the species that, when alone, 
maintains the greater equilibrium density will always displace i ts competitor 



(jAmarasekarel 120031 lO'Mallev et all . l2006al: lAllstadt et all . 12007^ 



Discrete, stochastic spatial models and their deterministic analogues have 
been employed, commonly and suc cessfully, to gain insight into collective beh av- 
ior of m ulti- "species" interactions ( Marro and Dickmanl 19991: Murravl 2003 ) in 
physi cs (korniss et al.l . ll995ill997() . chemistry ( Antal et al. . 19961 : Toroczkai et al 



1997; Zif f et al.l. 119861). and in the study of population dynamics (jEllner et al 



1998i: iMcKane and Newmanl . 12004 lO'Mallev et"al] . l2006bi l2009l) . We use both 
methods to investigate the combined effects of preemptive and interference com- 
petition on invasion. We organize our paper as follows. First, we present a 
discrete (individual-based), stochastic model where an invader and a resident 
species compete preemptively, and one or both species also employs interference. 
We let a species' propagation rate depend functionally on its level of interference; 
we consider convex, linear and concave trade-offs. We explore the model by an- 
alyzing invasibility criteria of both a mean-field approximation (homogeneous 
mixing) and a pair approximation. Then we apply results of the approximations 
to interpret simulations of the full spatial model. The Discussion compares our 
results to other spatial models incorporating allelopathy. Appendices collect 
much of the analytical detail. 

Our results find that interference by the invader increases its likelihood of 
successful invasion only when interactions are spatially structured. That is, the 
novel weapons advantage, in our model, appears only when we account for local 
clustering of invader individuals. Interference by the resident, the initially com- 
mon species, can inhibit invasion with and without spatial structure. Increases 
in background mortality rate (mortality before adjustment due to interference) 
diminish the competitiveness of the species relying more on interference. In- 
creased mortality reduces both local and global densities; that is, the frequency 
of empty sites increases. Consequently, the value of interference, relative to 
propagation, declines. Our model assumes that interference increases mortality 
of a preemptive competitor, but interference does not generate an alternative 
niche. Consequently, we should not anticipate co existence absent con tinuous 
introduction or strong effects of spatial clustering ( Allstadt et al. . 20091 ). 



2. A discrete, stochastic spatial model 

When plant species compete, inter actions regulating population growth gen- 



erally occur at the neighborhood scale ( Goldbera . ll987l:lUriarte et al. . 2004). In- 



terference competition, including allelopathy, ordinarily has a local spatial struc- 
ture. And, when propagule dispersal distance (inter-ramet distance in clonal 
species) is limited, plants compete preemptively for space at the neighborhood 
scale. Consequently, our model - which integrates preemptive and interference 
competition - assumes that competitive interactions a mong nearest neighbor s 
drive invasion and the dynamics of species' abundances (|0'Malley et al.l . l2006al ). 
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2.1. Model construction 

Two clonal plant species compete on an x Ly lattice with periodic bound- 
aries. Each lattice site represents the resources required to sustain a single in- 
dividual (a ramet) of either species. The local occupation number at site x 
is ?T.i(x) = 0, 1 with i = 1,2, referring to the resident and invader species, re- 
spectively. During a single simulated time unit, one Monte Carlo step per site 
[MCSS], Lr^Ly sites are chosen randomly for updating. 

An empty site may be occupied by species i through introduction from out- 
side the environment, or through local propagation. Introduction of species i 
at an open site occurs as a Poisson process with rate /3. Each species has the 
same introduction rate, to avert any effect of propagule pressure on the out- 
come of competition. Local propagation into an open site has rate ai77j(x), 
where ai is the individual- level propagation rate for species i, and r]i{x) = 
(l/(5)S]x'enn(x)'T'i(x') is the density of species i in the neighborhood around open 
site X. nn(x) is the set of nearest neighbors of site x, and S is the number of 
sites in that neighborhood (S — |nn(x)|). Since we equate colonization with 
propagation of new ramets, we let 5 = 4. 

Colonization can occur only at open sites. An occupied site opens through 
mortali ty. Individuals of each species suffer density-independent mortality at 
rate /z ( Cain et al.l . llQQSl) . An individual occupying site x experiences an in- 



creased mortality rate due to interference if nn{x) includes any heterospecifics. 
That is, an individual of species i at site x has total mortality rate ii+6jrij (x), i ^ 
J. 6j > represents interference by species j, and '7j(x) is the density of species 
j on the neighborhood around site x. 

Summarizing transition rules for an arbitrary site x, we have 

Q ;3+ au)l (x) ^ ^ /3+ a2r)2 (x) ^ ^ M+92»72(x) ^ ^ ^ 

where 0, 1, 2 indicates whether a site is open, resident-occupied, or invader- 
occupied, respectively. Table [T] defines the symbols and notation we use. We 
assume that interspecific competition drives the dynamics; i.e., each species per- 
sists absent competition. Therefore, we restrict attention to the (3 <^ ctdfi) < ai 
(i = 1,2) regime, where adu) is the critical propagation rate below which 
either species, in the other's absen c e, grow s too slowly to avoid extinction 
(lObornv et al.L l2005t lO'Mallev et"aLl . l2006al ). 



2.2. Life-history constraint 

The essential lesson of life-history theory is that increased allocation of re- 
sources to one tr ait advancing survival or rep roduction comes at the expense 



sources to one tr ait aavancmg survival or rep roduction comes at tne expense 
of another trait ( Bell and Koufopanoul . Il986l ). Applying this concept, we as- 



sume that any increase in a species' level of interference reduces that species' 
clonal propagation rate. The functional dependence (trade-off) has the form 
af + ^'i^ — C^, where C is a constant, equal to the maximal feasible rate of 



1^ 



ropagation, and R defines the shape of the trade-off (jCiraldeau and Caraco 



2000.) . If < i? < 1, the cost of increasing a (i.e., the decrease in 0) decreases 
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Table 1: Definitions of model variables and parameters 



Symbols 


Definitions 




Lx, Ly{— L) 


Lattice size 




X 


Location of lattice site 




ni(x) 


Occupation number for residents at site x; ni(x) = 


= 0,1 


n2(x) 


Occupation number for invaders at site x; n2(x) — 


0,1 


nii(x) 


Set of nearest neighbors around site x 




S 


size of neighborhood around site x ((5 = nn(x)|) 




f3 


Common introduction rate at empty sites 






Density of species i on nn(x) 






Individual rate of propagule production, species i 






Background mortality rate, both species 






Mortality of species i due to interference by species j 


Minimal propagation rate for persistence 




C 


Maximal propagation rate; C = 0.8 




R 


Curvature oi a,9 trade-off 




* 

Pi 


Equilibrium single-species density 





as a increases. If i? = 1, the cost of increasing a is constant, and if i? > 1, 
the cost of increasing a increases at greater a. More importantly, any increase 
in R increases competitiveness of a species that employs both preemption and 
interference; the constraint moves farther from the origin as R increases. If a 
species does not exercise interference, we let its propagation rate vary on {ac{p), 
0.8]. Therefore, we let C = 0.8. 

2.3. Novel weapons 

Under the novel weapons hypothesis, the invader exercises interference com- 
petition while the resident does not. That is, 6i = 0, and 



< 02 < ^C^—a,{f,f 



Invasive plants may interfere with native competitors through several mech- 
anisms. Some impact natives directly; others are mediated through a third 



speci es. Invaders can act as a disease reservoir (E ppinga et al 



2007 ). focus herbivory on native species (iDangremond et al 



2006; Borer et al 



201C)f ) , or p roduce 



allelopathic chemicals ( Prati and Bossdori , 20041 : Cipollini et al.l , 2008fl . Our 



analysis emphasizes the invader's trade-off between propagation and interference 
across a range of native species lacking interference. Under the novel weapons 
hypothesis, the resident species does not compete through interference, and we 
let its reproductive rate vary to model variation in the strength of preemption 
the invader encounters. 

2.4. Resistance zones 

To study resistance zones, we let the resident exercise interference competi- 



tion while the invader does not. That is, < 6*1 < yC^ — adfJ-) , and 62 = 
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A resident species exercising interference competition may strongly resist inva- 
sion since its initial density (by definition) will be relatively high. Given biotic 
resistance combining preemption and interference, we vary the invader's propa- 
gation rate to consider a range of matches between the invader and the abiotic 
environment invaded. 

2.5. Mutual interference 

Under the assumption of mutual interference competition, we have < 

61,62 < \Jc^ — ac{lj)^ . For this case we assume that each species' resource 
allocation between propagation and interference follows the same trade-off. 



3. Mean-field approximation 

The spatially-homogeneous mean-field approximation to our model ignores 
any effects of locally clustered growth. The mean-field model (MF) so off ers com- 
pariso n to results including inter act ions at the neighborhood scale : see WilsonI 



( 1998h . Pascual and Levin or lCuddington and YodzisI (|2onnl ) "or perspec- 



tive. 

pi and p2 represent the global densities of species 1 and species 2, respec- 
tively. Ignoring continuous introduction (letting j3 — Q) we have a mean-field 
dynamics: 

pi = ctipi{l ~ pi - P2) ~ pi{p + 62P2) (2) 
P2 = Ct2P2{^ - Pi - P2) - p2{p + 61P1) (3) 

6i represents the increased mortality of species j induced by species i\ we scale 
6i per unit density of species i. Species' persistence absent competition in the 
MF approximation requires only that each Ui > ac{p) — p. Letting 6i > 0, for 
i = 1, 2, we conduct a general stability analysis of the MF model's equilibria. 
Thereafter, we consider competitive superiority when the specific propagation- 
interference constraint introduced above applies. 

3.1. Equilibria and stability 

The dynamics has three boundary equilibria and one positive equilibrium. 
Mutual extinction, designated equilibrium El, cannot be stable since each > 
fi. At Equilibrium E2, species 1 competitively excludes species 2: 

E2: (^pi* = 1- =0^ (4) 

Species 2 excludes species 1 at equilibrium E3: 

E3: (pi* =Q,p2* ^1--^) (5) 
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Finally, a positive (internal) equilibrium, EA, can be expressed as: 



ai0i + 0262 + 9i02 ^ ' 

* _ 6>i(qi - /i) + - «2) 
" ai^i + 02^2 + 0162 ^ ' 



As noted above, we know that when competition is strictly preemptive (^i = 
02 = 0), the species cannot coexist under homogeneous mixing. 

Suppose that the resident excludes the invader. From Appendix A, local 
stability of invader extinction requires: 

+ <0 (8) 

ai ai 

Since /i < ai, 9i{-^ — 1) < 0. Therefore, any interference by the resident 
tends to stabilize invader extinction by increasing invader mortality. That is, 
interference by the common species can help prevent advance of the rare species, 
even if the rare species has the greater propagation rate. Clearly, if ai > 02 
and 9i > 0, the resident species repels the invader. Note that 62, the invader's 
level of interference competition, does not appear in Inequality dS]). Hence the 
invader cannot increase its growth rate when rare through interference; only 
increased propagation promotes mean-field invasion. 

To elaborate, the invader advances from rarity only if Inequality ([8|) is re- 
versed. From Appendix A, assuming 61 > 0, the rare species invades iff: 

a2 -ai > - 1) > (9) 

To invade successfully, the rare species must have the greater propagation rate, 
and must not suffer too much interference from the common species. 

Appendix A shows that if species 2 can invade the resident, then equilibrium 
E3 is locally stable. That is, if species 2 can advance when rare, it will grow to 
exclude species 1. 



3.2. Bistability 

Bistability requires that local-stability conditions for both single-species equi- 
librium nodes {E2 and i?3) hold simultaneously. Given bistability, initial condi- 
tions determine the outcome; at some point, the more abundant species "wins." 
From Appendix A, the common species (1 or 2) repels its rare competitor if 

9i(l -)<ai-a2 and 6l2(— - 1) > ai - 02 (10) 

Ai 

If ai > a2, the first expression must hold, since 6*1(1 — ^) < 0. The second 
inequality can hold simultaneously if 02 is relatively large, and the difference 
between propagation rates is not too large. If 02 > ai, the second expression 
must hold, since 62{— — 1) > 0, and symmetric conditions promoting bistability 
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are clear. Bistability can arise even if the species with the greater propagation 
rate does not exert interference. For example, if ai > a2, the system can be 
bistable when 9i = (provided, of course, 62 > 0). 

Figure [T] shows how the mean- field dynamics can flow in in the (pi, P2) phase 
space. One example plots competitive asymmetry. Species 2 advances when 
rare, and when common, excludes species 1. The other plots show bistability, 
one for a symmetric (identical species), and one for an asymmetric, domain of 
attraction. 

3.3. Coexistence? 

Coexistence (local stability of a positive equilibrium) requires that each 
species invade the other. Together, the conditions for mutual invasion indicate 
that coexistence requires (Appendix A): 



Since each > /z, each (1 — ^) < 0; the LHS of each of these inequalities 
must then be non-positive. But one (a^ — aj) must be positive, so that the two 
inequalities cannot be true simultaneously. Hence the MF does not permit mu- 
tual invasion, and so does not admit competitive coexistence, when the species 
interact through both preemption and interference. 

3.4. Propagation-interference constraint 

The preceding, general stability analyses did not treat interference as a func- 
tion of propagation rate. Invoking the life-history constraint generates a special 
case of the MF analysis. The condition for the invader's advance when rare. 
Inequality becomes: 



From the general MF analysis, any allocation to interference by the invader di- 
minishes its propagation rate (a2), and so decreases the likelihood of successful 
invasion. Any increase in R, relaxing the propagation-interference tradeoff, in- 
creases the range of feasible parameter combinations where the resident repels 
the invader (see fig. (H). That is, the opportunity for successful invasion de- 
clines as R increases in the MF model, since only the resident can benefit from 
interference under homogeneous mixing. 

Increasing background mortality (fi) increases the likelihood of successful 
invasion. Greater mortality, with di and the ai fixed, decreases the resident's 
density. Consequently, sites available for colonization increase in density, and 
the impact of interference on the invader's dynamics declines. Therefore, greater 
background mortality, which affects both species, decreases the range of feasible 
parameter combinations where the resident repels the invader (see fig. [2]). 

If the competitors mix homogeneously, interference can help the resident 
repel the invader, but cannot promote the rare species' invasion. Interference, 




(11) 




(12) 
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as a mechanism of biotic resistance, has a greater effect when a small level 
of interference does not cost too much in propagation, and when background 
mortality is relatively low (so that resident density is relatively high). 



4. Pair approximation 

Pair approximation (PA) incorporates corr elations of the occupation status 
of nearest neighboring sites into a dynamics ( Pickman . 19861 Matsuda et al 



[l987„; Bauch and Rand. 2000t | van Baalen . 2000,) . By incorporating a minimal 
spatial structure, PA predicts equilibr ium densities of individual-based model s 
more accurately than do MF models (l EUner et al. . 1998t Caraco et al. . 200 ll) . 



For our purposes, PA addresses consequences of local clustering, generated by 
dispersal limitation, for invasion. Appendix B presents details. 

The PA tracks global densities Pi{t), and conditional densities qj\i{t) rep- 
resenting the likelihood that a site is in state ?', g iven that a neighboring site 



resenting tne liitelinooa tnat a site is m state ?, g iven tnat a neigiiDormg site 
in state i (jSato et all . l2000t [o^allev et all l2006al) . The three global densities 



{po{t), pi{t), p2{t)) imply 9 local densities {qj\i)- However, the PA's dimension 
is limited by simple constraints: 

Po + Pi + P2 = 1 
9oK + 9i|j + 92K = 1 (13) 

1i\]P3 = <lj\tPii 

where i,j = 0,l,2;i ^ j; wc suppress time dependence for simplicity. These 
constraints leave only 5 of the 12 total variables are independent. Since species 
2 is the invader, we track pi, P2, '?2|2j 9112; ^ili- Note that the dynamics of qi\2, 
the conditional density of a resident given an invader at a neighboring site, can 
depend on a novel- weapons effect when 6*2 > 0. 

The dynamics of the global densities, like the mean field approximation, 
account introduction to empty sites, local propagation, background mortality, 
and death due to interference competition. For the resident, 

Pi = /3po + Piai^oli - PiM - ^'2'72|iPi (14) 

Compared to the mean-field model, both local propagation and mortality due 
to interference depend on local, rather than global, densities. Hence the PA 
models both site preemption and interference as effects of spatially clustered 
growth. 

To model the dynamics of a local density, we first write the dynamics of a 
doublet Pa, a global density, where pn — Piqi\i. The PA's doublet dynamics 
{dpii/dt; see Appendix B) introduces the triplets g2|02 ^nd 91122- Ordinary 
PA assumes that neighbors of neighboring sites are weakly correlated, and lets 
92|02 = 92|o and qi\22 = qi\2- The approxima tion closes the syst em of equations, 
permitting us to write an invasion criterion ( Iwasa et al.l , 119981) . 



For the invasion analysis, assume that an introduction of species 2 has oc- 
curred. Invasion either succeeds or fails before the next introduction event 
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occurs (since /3 <^ fi < ai). Successful invasion requires that the invader have 
a positive growth rate when rare; p2 > as p2 ^ 0. From Appendix B, this 
condition yields the invasion criterion: 

l-92|2-9i|2> ■ (15) 

The left side represents the frequency of open sites neighboring an invader, 
go|2- The right side represents the ratio of death rate to birth rate for the 
invader. Successful invasion, then, requires that the invader colonize a neigh- 
boring, empty site before it dies. The death rate sums background mortality 
and averaged mortality from the resistance zone about a resident neighboring 
the invader. The seeming simplicity of the invasion criterion masks an impor- 
tant biological difference between the MF criterion and the neighborhood-scale 
condition for invasion. In the PA condition for successful invasion, the local 
density gi|2, hence the local density of open sites go|2j depends on 02, invader 
interference (see Appendix B). That is, the neighborhood-scale invasion criterion 
reveals a role for novel weapons, contrasting to the MF result. Consequently, an 
interfering invader might invade a resident species despite the resident having 
the greater propagation rate. The right side of the criterion, of course, include 
the resistance-zone effect of interference by the resident. 

To compare the PA invasion criterion to both the MF and simulation of the 
full spatial model (see below), we evaluated Inequality (|15l) numerically. We 
constructed pairwise invasion plots for the three cases defined above: only the 
invader interferes (novel weapons) , only the resident interferes (resistance zones) 
and both species exercise interference competition. Whenever a species was an 
interference competitor, we invoked the propagation-interference constraint. 

Figure 13] shows pairwise invasion results when the invader, and not the resi- 
dent, is the interference competitor. That is, the invader obeys the propagation- 
interference constraint, and 9i = independently of the resident's propagation 
rate; we so isolate the novel weapons effect. When interactions are local, in- 
terference competition can promote an invader's growth when rare. When the 
propagation cost of interference is smaller (R is larger) and ^ is small, invaders 
with a propagation rate much less than the resident's rate succeed. Invader in- 
terference proves advantageous since it opens sites neighboring invaders. Hence 
interference makes space available where the invader's local density may match 
or exceed the resident's local density - despite the resident's greater global den- 
sity. Hence, novel weapons, in our model, enhances invasion when individuals 
compete at the neighborhood scale, but has no effect under homogeneous mix- 
ing. 

Figure [3] shows pairwise invasion results when the resident, and not the in- 
vader, is the interference competitor. That is, the resident obeys the propagation- 
interference constraint, and 62 = independently of the invader's propagation 
rate; we isolate the resistance zone effect. A resident can repel an invader with a 
much higher propagation rate than its own through local interference. Increas- 
ing R now reduces the range of parameter combinations permitting invasion, and 
increasing fi promotes invader success - results symmetrically opposite to those 
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for novel weapons. When only the resident can interfere, a greater value for R 
reduces the pleiotropic cost of interference. Increasing ^ opens more space for 
the invader; greater background mortality renders a given level of interference 
less effective as a competitive mechanism. 

Figure [S] shows pairwise invasion results when both species exercise interfer- 
ence, and so both are constrained by the propagation-interference tradeoff. Not 
surprisingly, the resident species, because of its initial high density, can employ 
interference to repel the invader over a much greater range of parameters than 
when both species lack interference (where greater propagation always wins). 
The utility of interference, from the resident species' perspective, declines at 
small R, and at greater fi; the invasion plot for (i? = 0.5, /x = 0.2) matches the 
case where competition is strictly preemptive. In most of the plots substantial 
regions of bistability appear. A species with intermediate levels of both prop- 
agation and interference can repel a broad range of propagation-interference 
combinations. And, when rare, the intermediate combination is repelled by 
species within that broad range, as long as i? > 1. 

Mutual interference does not permit coexistence under PA (fig. [5]). In gen- 
eral, the invasion plot resembles the case where only the resident can interfere. 
But for sufficiently large R and small enough fj,, a resident that invests solely 
in propagation, and not in interference, can be invaded by a species that mixes 
propagation and interference (rightmost columns of invasion plots). The resi- 
dent, with an initially high density and no interference capacity, is vulnerable 
to a clustered invader that can locally open sites via interference competition. 



5. Simulation of individual-based model 



We simulated the individual-based spatial model with ^ Ly = L = 256, 
and (3 = 0.001. We initiated simulations with the resident occupying each site. 
The resident's density was allowed to decline to its single species equilibrium 
density, pj, without any introduction. Then, at time t = 0, introduction of 
both invader and resident individuals began. We tracked the global densities of 
each species pi (t) , and recorded a successful invasion when the invader reduced 
the global density of the resident species to pl/2 within 20,000 MCSS time 
steps. These simulations envision growth of initially small invader clusters at 
the initiation of ecological i nvasion ( Korniss and Caracol 2005; O'Malley et al 
2006at lAllstadt et al. 1 . 120071 ) . 



We can attribute some differences between results of the individual-based 
spatial model and its pair approximation to the difference in degree of spatial 
clustering. PA, by construction, truncates correlation length at a nearest neigh- 
bors. The simulation model permits development of greater correlation lengths 
(invader clusters can grow large from a single site). Cluster size impacts the 
relative frequency at whi ch individuals experien ce local self- regulation versus 



interspecific competition ( O'Mallev et al. . 2010t ). These frequencies, in turn 



affect invasion success and global population densities. 

We simulated invasion from rare introduction for the three cases examined 
above: novel weapons, resistance zones, and the case where both species interfere 
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competitively. Figure [S] shows pairwise invasion results when the invader, and 
not the resident, interferes competitively in the individual-based spatial model. 
The invader obeys the propagation-interference constraint, and 9i — inde- 
pendently of the resident's propagation rate, isolating the novel weapons effect. 
Invader interference promotes successful invasion more frequently in simulation 
than under pair approximation. That is, we note a novel-weapons effect (absent 
in the MF) . Invader clustering increases the advantage of local interference com- 
petition against a resident lacking interference. Even when the propagation-cost 
of interference is steep {R < 1) an invader may advance and exclude a resident 
with a greater propagation rate. When interference is less costly {R > 1) most 
invaders succeed, except those with a very low propagation rate. PA, then gives 
a reasonable, but understated, prediction of the impact of novel weapons in the 
detailed spatial model. 

Figure [7] shows pairwise invasion results when the resident, and not the 
invader, interferes competitively in the individual- based model. The resident 
obeys the propagation- interference constraint, and 02 = independently of the 
invader's propagation rate, isolating the resistance zone effect. Interference 
allows a resident to repel invasion by rare species with much greater rates of 
local propagation. For the resistance-zone effect, PA predicts the results of the 
individual-based model very accurately across all parameter combinations we 
simulated. 

The simulation model's pairwise invasion results for the case where both 
species exercise interference competition appear in Figure [H When the cost of 
interference is high (i? — 0.5) the results are very close to those for preemptive 
competition only. That is, the species with the greater rate of propagation wins 
in most cases. Compared to PA, invader clustering in the simulation model pro- 
duces more cases of successful invasion at low background mortality. At lesser 
costs of interference {R = 1,2) residents mixing intermediate levels of propaga- 
tion and interference repel most invading species, particularly when background 
mortality is not too large. Residents that invest either too much in interference 
(low ai) or too little are susceptible to invasion by species with a more balanced 
mix of propagation and interference. 

In a separate exercise, we released both species from the propagation-interference 
constraint and searched the parameter space for competitive coexistence. We set 
/i = 0.1, /? = 0.001, and let i = 64 to save computation time. The simulations 
ran for 100,000 MCSS, and the threshold for coexistence was pi > 0.05; i = 1,2, 
well above the expected density due to introduction alone. The search yielded 
no evidence of competitive coexistence, as our analytical models predict. 



6. Discussion 



In the novel weapons scenario an invader gains competitive superiority over 
a resident through interference. A pplications focus on allelopathic interference 
( Callaway and Ascheghouj . 2000t ). If species mix homogeneously, we find that a 



rare i nvader gains no advantage through interference competition ([Case and Gilpin 
1974 . But if competitive interactions occur at the local scale, an invader can 
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employ interference to advance against a resident species with a greater rate of 
local propagation (figs. [3]and[6]). Invaders are rare globally but locally clustered 
in our spatial models. Invaders can use interference to reduce resident density 
at the perimeter of these clusters, wh ere local invader dens i ty is s ufficiently high 
to compete effectively for open sites (jKorniss and Caracoll2005h . 

In the resistance zone scenario, the resident holds competitive superiority 
over an invader through interference. A resident's competitive interference has 
a strong effect under homogeneous mixing; few invaders advance successfully 
unless interference is costly and background mortality is high (fig. Under 
pair approximation, interference by the resident restricts invasion in two cases: 
when the resident, but not the invader, interferes, and when both species exercise 
interference competition. The resistance zone effect also appears in simulation 
of the individual-based model. But when both species can interfere, successful 
invasion is far more common in the simulation. The difference between the full 
spatial model and its pair approximation lies in the impact of invader clusters 
larger than neighborhood size (hence, longer correlation distances) generated 
by the individual-base model. Clustering increases an invader's benefit from 
interference against an interfering resident. 

When both species exhibit interference, the pairwise invasion plots indicate 
regions of bistability (if only one species interferes, one cannot meaningfully 
reflect invasion plots about the diagonal). We found large regions of bistability 
under homogeneous mi xing. Bistability declined as spatial structure increased 
(jChao and Levinl . Il98l[) , and nearly disappeared in the simulation model (fig. 
[8]). Given a sufficiently long time, we would anticipate that the spatial process 
on a large, but finite, lattice would result in competitive exclusion. The time 
required, however, may exceed population-dynamic scales. 

In a three-species spatial model, where a non-int erfering species interacted 
with both a weak interferer and a strong interferer, iDurrett and Levinl ()l997l) 
found cyclic coexistence. Our models for two-species competition, by construc- 
tion, do not admit coexistence. In fact, when we assumed that both species 
exercise interference competition, and when increased interference reduces prop- 
agation, certain parameter sets suggested that a single combination of propaga- 
tion and interference might repel all other feasible combinations. Suppose we 
start the resident's propagation rate, ai, at minimal (maximal) levels. Then 
larger (smaller) propagation rates successively invade and exclude the resident 
until pairwise competition leaves a resident species that can exclude all others. 
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Appendix A. Mean-field stability 

The Jacobian for the spatially homogeneous MF dynamics is: 



J = 



ai - /X - 2aipi* - p2*(ai + O2) -pi*{ai + 62) 

-P2*{a2 + di) a2 - M - 2a2P2* - Pi*(a2 + di) 



Evaluated at El, mutual extinction, the Jacobian becomes: 
J{E1) = 



ai — p, 
a2 - p 



(A.l) 
(A.2) 



Since ai > p for i ~ 1, 2, by assumption, mutual extinction is unstable. 
At E2 the resident excludes the invader. The associated Jacobian yields: 



J{E2) = 
The two eigenvalues are 



a2-M-(l-^)(a2 + 0i) 



Ai(i?2)=/i(^-l) + 0i(-^-l) 
ai a\ 

\2{E2) =p-ai<0 



(A.3) 



(A.4) 
(A5) 



Local stability requires Xi{E2) < 0. 

For species 2 to advance when rare, Ai(i?2) must be positive (implying that 
extinction of species 2 is unstable). Note that if 02 > ai , then — 1) > 0, a 
condition promoting invasion by species 2 when species 1 rests its single-species 
equilibrium E2. When 9i > 0, species 2 invades iff: 



a2 — ai > 9i{ 



Oi 



1) > 



(A.6) 



Next, consider E3, where species 2 competitively excludes species 1. Evalu- 
ating the Jacobian, we have: 



J{E3) 



- p - (l - £){a, + 62) 
-(l-£)(«2 + ei) a2-At-2a2(l 



02 ' 



(A.7) 



The symmetry between J{E2) and J{E3) extends to the eigenvalues. From 
Eg IA.7[ we obtain: 



XiiE3)^pi 



^-l) + ^2(^-l) 
a2 a2 



X2iE3)=p-a2 <0 
Species 2 cannot exclude species 1 ijf 

ai - a2 > 02{. 1) > 



(A.8) 
(A9) 

(AlO) 
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Clearly, reversing subscripts in Expression ([SJ yields Expression (|A.10[) . More 
importantly, we see that if Inequality ([9]) holds, then Ai(£'3) < 0. Hence, if 
species 2 can invade species 1, i?3 is the only stable equilibrium, and species 2 
advances to exclude species 1 competitively. 

Now consider bistability. When common, species 1 repels invasion by species 
2 if Xi{E2) < 0. When species 2 is common, it repels invasion by species 1 if 
Ai(i?3) < 0. Both conditions hold, and the dynamics is bistable, iff: 

ei{l -)<ai-a2 and 6l2(— - 1) > ai - (A.ll) 

From our analysis of the single-species equilibria, the two species can coexist iff 
Xi{E2) > and Ai(i?3) > 0. The first inequality implies that when species 1 
rests at its single-species equilibrium, species 2 can advance from rarity. The 
second inequality reverses roles of common and rare. In the text we show that 
the mean-field dynamics does not admit competitive coexistence. 



Appendix B. Pair-correlation dynamics 

We deve l op ou r pair approximation (PA) by modifying methods described by 
Iwasa et al.l ( 19981 ). We write a dynamics for five state variables: pi, P2,q2\2T'li\2TQi\i- 



The first two are global densities, and the other three are local densities. In- 
voking constraints listed in the text. Expression (|14l) . we express the remaining 
PA variables in terms of the five state variables. In particular: 

Po = l-Pi-P2 1i\a = -. (1 - - 9i|2— j 

i — Pl — P2 Pi 

92|0 = (-; ) (1 - 92|2 - 9112) 90|2 1 - 9212 ~ qi\2 

^ — Pi — P2 

-, P2 P2 

9o|i = 1 - - 9i|2— 92|i=9i|2 — 

(B.l) 

We omit qo\Q, also easily calculated, since it does not appear in the analysis. 

To begin, we rewrite the dynamics of the resident's global density, Eq. ([Ti)) . 
in terms of our five state variables. Doing the same for the invader's global 
density yields: 

Pi = ^(1 - Pi - P2) + Pi [ai(l - gill - gi|2 — ) - P - ^'2'7i|2— ] (B.2) 

Pi Pi 

P2 = f3{l - pi - P2) + P2 ["2(1 - q2\2 - qi\2) - P - ^'19112] (B.3) 

To write the dynamics of the three conditional densities, we first require the 

Pi 



dynamics of a corresponding doublet. By definition, qj\i = where pij is the 
unordered doublet density. The doublet pij has dynamics: 



qJ\^ = -P.J - —P^ (B.4) 
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Applying the recipe to q2\2, we have: 

P22 = 2po2^ + 2po2Q!2 7 ^ ^92102 



2p2 



p 



-6*1511 



22 



(B.5) 



The first two terms represent generation of new invader pairs through intro- 
duction and birth, and the third term represents the loss of invader pairs due 
to background mortahty and interference competition. This equation intro- 
duces triplets 92102 a-nd qi\22 into the dynamics. Ordinary PA takes neighbors of 
neighboring sites as weakly correlated, and so we assume that 92102 = 92 |o a-nd 
9i|22 = 9i|2- The resulting closure of the equations allows the analysis without 
including the 27 types of triplets, or any higher order spatial correlations. 

Using the closure assumption and converting terms with equation (IB.ip . 
equation (|B.5p becomes 



P22 = 2p2(l - 92|2 - ?1|2) /3 + a2(^ + T7Y~r— - 9212 - ?1|2)) 



'2p292|2 



(5(1 -pi- P2) 

5-1 
P H ^6'i9i|2 



(B.6) 



Substituting Eqq. and (|B.6|) into Eq. (|B.4|) yields: 



92|2 = 2(1 - 92|2 - 9112) 



a2 



1 



P2{5-1) 



(1 - 92|2 - 



-92|2 



P2 



5 ^ 1 - Pi - P2 

(1 - pi - P2) + 02(1 - 92|2 - 9112) - P 

5-1 



9112)) 

6*19112 



-292|2 



p 



7l9l|2 



(B.7) 



Similarly, 



?i|i = 2(1 - 9iM - 9i 



'Pi' 



-9i|i 



^ + ^(1 + -; (1 - gi|i - 9i|2 — )) 

1 - Pi - P2 Pi 



(1 - pi - P2) + ai(l - 9i|i - 9i|2 — ) - P~ 6*29112 



-291 



P- 



Pi 

Ifl P2 

-"29112 

Pi 



Pi 



(B.E 



P12 is slightly more complicated. There are two ways to make a (1, 2) pair, and 
each member of a (1, 2) pair interferes with the other. Proceeding: 



P12 = P2(l - 92|2 - ?1|2) 



-Pi(l - gill - ?i|2 — ) 
Pi 



aipi{5-l) p2, 

Yf. r(l - gi|i - gi|2 — , 

0(1 - Pi - P2) Pi 



/3 



a2P2{5 - 1) 
(5(1 - pi - P2; 



(1 - 92|2 - 9112) 



-P29l|2 



2p+^{l + {5- 1)9112) + -f (1 + (5 - l)9i|2-; 
d ' Pi 



(B.9) 
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Then, following substitution: 



9112 = (1 - 9212 - qi\2) P + T (1 ^ ^ Il\2 — J 

' L d(l - Pi - Pa) Pi 



+Pi(l - 

-gi|2 



9l|2 



Pi' 



/3 



a2('5-l) 



-9i|2 



P2 (5(l-pi-p2) 

/3(1-Pi-P2) 



(1 - 92|2 - gi|2y 
92 



A^ + f (l + (5-l)9i|2) + f (l + (<5-lkl|2^) 



P2 



+ 02(1 - 9212 - 9112) - 911 



2^1 



(B.IO) 



Equations (jB.2p . (B.3), (jB.7p . (jB.8|) . and (|B.10p constitute the pair- approximation 
dynamics. 

To analyze invasion, we introduce species 2 at near-zero density. Invasion 
succeeds or fails before the next introduction event occurs (since /3 <C /x^ < Qi). 
Taking /3 = 0, the PA dynamics becomes: 



pi = pi [ai(l - qi\i - gi|2 — ) - ^J'- 6*29112—] 

pi pi 

P2 P2 ["2(1 - 92|2 - 9112) - M - 6*19112] 



P2l 



92|2 = 202(1 - 92|2 - 9112) 

-9212 



9i|i = 2ai(l - qi\i - 9i|2— ) 
' ' pi 



ri , P2(^-i) 
(5 5(1 - pi - P2) 

l|2)+M + 

Pi(5-1) 



(1 - 921: 



9i|2; 



"2(1 - 92|2 - 9112) + M + 6*19112 (^^^-^-^ - 1) 



5 



1 , PU" - ^ /, P2, 

T + TTj r(l - 9i|i - 9i|2 — ; 

d(l - pi - P2) Pi 



-9i|i 



n P^\^ P2, 2(<5-1) 

ai(l - 9i|i - 9i|2 — ) + M + 6'29i|2 — ( j 

Pi ' pi ' d 



1 



9i|2 = (1 - 92|2 - 9i|2)(l - 9i|i - 9i|2— )(t7^^^^^ ^^)(ai + "2) 

' pi 'd(l - pi - P2)' 

-9i|2 A* + "2(1 -92|2 -9112) 



(B.ll) 
(B.12) 
(B.13) 

(B.14) 
(B.15) 



-9i|2 



The next step of the invasion analysis addresses the frequency of open sites when 
the invader is rare. Since the invader is rare, species 2 has no effect on either the 
resident's equilibrium global density (pi) or the equilibrium frequency of paired 
residents (9i|i). Setting p2 = in Eqq. (|B.lip and (|B.14p yields: 



and 



Pi = 



9i|i 



ai 

(5- 1 - ^ 

ai 



ai 



(B.16) 
(B.17) 
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Now we use these results to find the other conditional probabilities. Let x = q2\2, 
y = 9i|2, and w = go|2 = I - x - y. Then: 



1 



2(6-1) 



11 - )[ax +a2)w + y 



1 



(5 

(5-1 



a2wx — fj,x 
01 6 



(B.18) 



Solving for the equilibria, we have: 



y 



a2Sw + ^iS + y*6i{5 - 2) 

w{ai + a2)((5 - 1 - 5-^) 
(5 [/i + a2w\ + 6*1 - y*ei + 6*2 



(B.19) 
(B.20) 



Given the interdependence of x* and y*, we evaluated {p2)p2=o > numerically. 
Using Eq. (jB.12|) . we obtain the invasion criterion, Eq. ([T5|) in the text. The 
importance of deriving x* and y* lies in demonstrating that at the neighborhood 
scale, both species' level of interference [Oi and 82) affects the likelihood that 
species 2 can invade the resident species 1. That is, the invasion criterion defined 
by neighborhood-scale correlations depends not only on the resident's level of 
interference 9i (resistance zone), but on the invader's level of interference 62 
(novel weapons) . 
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Figure 1: Flow in the mean-field dynamics over the (pi,P2) phase space, fj, = 0.25 in each 
plot. Top: Species 2 excludes species 1 independently of initial condition; Expression (|9]l 
holds, ai = 0.5, 9i = 0.2, and 02 = 0.75. Middle: Condition for bistability. Expression HIOI I. 
holds, ai = 0.5 = 02 = 0.6, 61 = 62 = 0.1. Bottom: Bistable dynamics; Qi = 0.5 = a2 = 0.6, 
61 = 0.2, 82 = 0.1. Species 2 has the greater domain of attraction. 
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Figure 2: Pairwise invasion plots generated by the mean-field approximation. Each subplot's 
ordinate is the resident's propagation rate, ai. The resident's level of interference is given 
by the life-history constraint, -|- 6^ = C^; The MF invasion criterion is insensitive to 
interference by the invader. Black: parameter combinations where invader's growth is positive 
when rare; invasion succeeds. White: resident repels invader; invasion fails. Separation of 
these equilibrium phases follows from Inequality II12I I. Note that the diagonal (solid line) 
separates invader success from failure when each species competes through site preemption 
only. Rows, top to bottom, show that invasion increases as background mortality increases. 
Columns, left to right, show that invader's success declines as resident propagation increases 
for given level of interference. If any of the plots is reflected along the diagonal, parameter 
combinations where "white falls on white" fulfill criteria for bistability. Dotted lines denote 
cec{fi) for each species, and neither species is viable in the hatched region. 
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Figure 3: Novel weapons under pair approximation. Pairwise invasion jjlots generated by 
the PA invasion criterion; only the invader exercises interference. Black: parameter combina- 
tions where invader's growth is positive when rare; invasion succeeds. White: resident repels 
invader; invasion fails. 
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Figure 4: Resistance zones under pair approximation. Pairwise invasion plots generated by 
the PA invasion criterion; only the resident exercises interference. Black: parameter combina- 
tions where invader's grovjfth is positive when rare; invasion succeeds. White: resident repels 
invader; invasion fails. 
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Figure 5: Pairwise invasion plots generated by the PA invasion criterion; both species exercise 
interference. Bla<;k: parameter combinations where invader's growth is positive when rare; 
invasion succeeds. White: resident repels invader; invasion fails. 
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Figure 6: Novel weapons in the simulation model of invader clusters. Each point represents 
combined results of 20 simulations. Bla<;k: parameter combinations where invader's growth is 
positive when rare; invasion succeeds. White: resident repels invader; invasion fails. Stochas- 
ticity led to mixed results for a few parameter combinations, which axe shaded gray according 
to the proportion of successful invasion. 
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Figure 7: Resistance zones in the simulation model of invader clusters. Each point represents 
combined results of 20 simulations. Black: parameter combinations where invader's growth is 
positive wrhcn rare; invasion succeeds. White: resident repels invader; invasion fails. Gray: as 
described above. 
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Figure 8: Simulation model: invader clusters when both species interfere. Each point rep- 
resents combined results of 20 simulations. Black: parameter combinations where invader's 
growth is positive when rare; invasion succeeds. White: resident repels invader; invasion fails. 
Gray: as described above. 
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